Relaxation dynamics of a linear molecule in a random static medium: 

A scaling analysis 
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We present extensive molecular dynamics simulations of the motion of a single linear rigid molecule 
in a two-dimensional random array of fixed overlapping disk-like obstacles. The diffusion constant 
for the center of mass translation, Dcm, and for rotation, Dr, are calculated for a wide range of 
the molecular length, L, and the density of obstacles, p. The obtained results follow a master curve 
Dp M ~ (L 2 p)~ v with an exponent fi — —3/4 and 1/4 for Dr and Dcm respectively, that can be 
deduced from simple scaling and kinematic arguments. The non-trivial positive exponent v shows 
an abrupt crossover at L? p — £i. For Dcm we find a second crossover at L 2 p — (2- The values of 
£1 and (^2 correspond to the average minor and major axis of the elliptic holes that characterize the 
random configuration of the obstacles. A violation of the Stokes-Einstein-Debye relation is observed 
for L 2 p > £1, in analogy with the phenomenon of enhanced translational diffusion observed in 
supercooled liquids close to the glass transition temperature. 
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I. INTRODUCTION 

Since it was initially introduced by Lorentz as a model 
for the electrical conductivity in metals [1], the problem 
of the Lorentz gas has given rise to a substantial theo- 
retical effort aimed to understand its properties [2-10]. 
In this model, a single classical particle moves through a 
disordered array of static obstacles. It can thus be used 
as a simplified picture of the motion of a light atom in a 
disordered environment of heavy particles having a much 
slower dynamics. It has been shown that, for high densi- 
ties of obstacles, this model presents many of the features 
of the dynamics of supercooled liquids or dense colloidal 
systems, such as a transition from an ergodic phase of 
non-zero diffusivity to a non-ergodic phase with diffusiv- 
ity zero [2-6,8]. Despite of the absence of a dynamics of 
the host medium, theoretical approaches to the Lorentz 
gas problem are highly non-trivial. In particular, dif- 
fusion constants and correlation functions are found to 
be non-analytical functions of the density of obstacles 
[2,4-6,9,10]. 

Diffusing particles and obstacles are generally mod- 
eled as disks, spheres or hyperspheres in two, three, and 
11 > 3 dimensions, respectively. In this paper we present 
an investigation by means of molecular dynamics simu- 
lations of the relaxation dynamics of a generalization of 
the Lorentz gas, in that the diffusing particle is given 
by a rigid rod. This model can thus be seen as a sim- 
plified picture of the dynamics of a linear molecule in 
a porous medium or in a colloidal suspension of heavy 
particles. Furthermore this system will allow to gain 
insight into the relaxation dynamics of a nonspherical 
probe molecule immersed in a simple liquid, an exper- 
imental technique that is used, e.g., in photo- bleaching 
experiments to study the dynamics of glass-forming liq- 
uids [11,12], or the investigation of the coupling between 
the rotational and translational degrees of freedom in 



such liquids [13,14]. 

For simplicity, the present simulations have been done 
in two dimensions, though we expect that physical ar- 
guments similar to those presented here can be used for 
the analysis of the results in the more general three- or 
n-dimcnsional case. 

The paper is organized as follows: The details of the 
investigated model and the simulation are given in Sec- 
tion II. Qualitative information about the trajectories of 
the rod is presented in Section III. Section IV presents 
results for the rotational and center-of-mass translational 
diffusion constants of the rod. An interpretation of the 
obtained results in terms of characteristic lengths of the 
host medium is given in Section V and the conclusions 
are presented in Section VI. 



II. MODEL AND DETAILS OF THE 
SIMULATION 

The generalization of the Lorentz gas we consider is 
a rigid rod diffusing in a two-dimensional array of static 
disk-like obstacles. The rod, of total mass M, consists 
of a string of N linearly aligned beads of equal mass 
to = M/N with a bond length 2a. The rod length 
is therefore given by L — (2N — l)a. Although from 
the point of view of a theoretical calculation the use of 
hard-core like interactions is preferable, it is simpler to 
do a simulation with a continuous potential. Thus, we 
modeled the interaction between single beads and sin- 
gle obstacles by a truncated and shifted soft-disk po- 



tential, V(r) = e[{(j/r) 12 — (a/r c ) 12 ] for r < r c , and 
zero for r > r c . As cutoff radius we took a value of 
r c = 2.5(7. The obstacles were placed randomly in a 
square box of length ^box, which was used for periodic 
boundary conditions. Obstacles were allowed to overlap, 
i.e., no minimal distance was imposed between the cen- 
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tcrs of neighbouring disks. The density of obstacles is 
defined as p — rc bs/^box> with n Q b s the number of obsta- 
cles. Subsequently a rod was randomly placed in the box 
and was equilibrated at a temperature T = e/fce for a 
time longer than the relaxation time for the correspond- 
ing p and L. (Note that this equilibration was necessary 
since interactions arc not hard-core.) In the following, 
space and time will be measured in the reduced units a 
and (cr 2 m/e) 1 / 2 , respectively. 

Once the rod was equilibrated, a production run was 
performed in the microcanonical ensemble. The equa- 
tions of motion of the rod were formulated in terms of 
the coordinates of its center of mass and its orientation. 
Due to the linear geometry of the rod, the constant length 
of the bonds and equal mass of the beads, the total force 
and the torque of the constraining forces on each single 
bead is zero, and thus only the total force and the torque 
due to the external forces of the obstacles on the beads 
are needed in the equations of motion. The latter were 
integrated by using the velocity form of the Vcrlct algo- 
rithm. We used a time step between 0.01 and 0.002, re- 
spectively for the smallest and the largest investigated L 
and p. Production runs covered typically about 10 6 time 
units. For each configuration of obstacles we simulated 
100 independent rods and in order to take the average 
over the frozen disorder we investigated typically 6-10 
different realizations of obstacles. Thus, for each density 
and rod length we simulated between 600-1000 indepen- 
dent rods. The number of obstacles was typically 50000. 

III. TRAJECTORIES 

Fig. 1 shows for different rod length typical trajecto- 
ries of the center of mass of the rod over a time of 10 6 
time units and the density p = 10~ 2 . For short rods, 
L = 3, the trajectory looks like the one expected for a 
diffusive Brownian particle. With increasing L the na- 
ture of the motion changes not only quantitatively but 
also qualitatively. As it is recognized from the figure for 
the case L — 79, for long rods the trajectory is composed 
of vertices in which the rod rattles back and forth and 
long segments in which moves in a quasi-linear way. The 
reason for this change in the dynamics is that, for large 
L, the rods can be temporarily trapped by the neighbor- 
ing obstacles, an effect that is well-known in the case of 
the original Lorentz gas for L = 1 and large p. For this 
latter system an increase of p leads to a strong increase 
of the trapping time and hence to a slowing down of the 
dynamics [2,3], in analogy to the mechanism that leads 
to the dramatic slowing down of the relaxation dynamics 
in simple glass-forming liquids upon cooling. This trap- 
ping mechanism is of course also present in the case of 
the rods. However, while for the original Lorentz gas, 
i.e., for disks, the diffusion dynamics is isotropic also on 
small length scales, the translational relaxation dynam- 
ics of rods in the direction perpendicular to their axis 



is strongly hindered by the steric hindrance induced by 
the neighboring obstacles, leading to the observed quasi- 
linear paths, in which the velocity of the rod is approxi- 
mately parallel to its axis and thus allowing for a quick 
propagation. Although these events are relatively rare, 
they allow the rod to move rapidly through the obstacles 
and hence they are important for its relaxation dynamics. 
Since the time it takes for the rod to have a noticeable 
velocity component that is (essentially) parallel to its ori- 
entation and the length of the subsequent quasi-ballistic 
flight is a random variable (due to the strong disorder of 
the obstacles), the resulting motion is very erratic and 
also strongly heterogeneous (i.e. it depends strongly on 
the initial position and orientation of the rod). Hence it 
is no surprise that the resulting average dynamics is also 
very heterogeneous which can be seen, e.g., by the pres- 
ence of strongly non-Gaussian effects in the relaxation 
dynamics at large length scales [15]. 

Fig. 1 also shows that the typical spatial extension 
of the trajectory within a given time depends strongly 
on L which implies that also the translational diffu- 
sion constant for the center of mass shows a strong 
L— dependence, in agreement with the results that we 
will discuss below. Hence it might be suspected that the 
change in the nature of the trajectories upon a change 
in L, see Fig. 1, is just due to the decrease in Dcm- 
That this is, however, not the case, can be seen in Fig. 2, 
where we show two typical trajectories for two rods that 
have a very different length but, by having made an ap- 
propriate choice of the corresponding densities p, a very 
similar translational diffusion constant for the center of 
mass. As can be seen, the nature of the trajectories does 
indeed depend strongly on p and L, even if the diffusion 
constant is the same. 



IV. DIFFUSION CONSTANTS AND SCALING 
LAWS 

We now investigate the L— and p— dependence of the 
diffusion constants for rotation, Dr, and for transla- 
tion of the center of mass, Dqm- These quantities 
can be calculated as Dr = ((A(f>) 2 } /2t s i m and Dqm = 
((Ar) 2 )/4t sim , respectively, where ((A</>) 2 ) and ((Ar) 2 ) 
are the mean square angular and center of mass transla- 
tional displacements at t = t s ; m , the time at the end of 
the simulation. 

We first consider Dr and we assume that its depen- 
dence on L and p is given by a product of the form 

D R = k R L- a p- b , (1) 

with k R a constant prefactor. By using scaling and 
kinematic arguments we show in the following that this 
Ansatz does indeed gives a good description of the data 
in that we can generate a master curve for D R (L, p). To 
this aim we note that in the limit of infinitely thin rods 
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and point obstacles, i.e., in the limit that finite size ef- 
fects are absent, the geometry of the system is invariant 
if it is scaled by an (arbitrary) factor /. Thus, since the 
length of the rod and the density of obstacles for the 
scaled system will be given respectively by fL and f~ 2 p, 
it is clear that the relevant quantity for the geometry of 
the dynamics is the dimensionless variable L 2 p, since it 
is invariant under scaling. It is therefore useful to rewrite 
the Ansatz for D R in Eq. (1) as: 

At = k RP a /(L 2 pf . (2) 

The corresponding rotational diffusion constant for a 
given system, D R , and that of the rescaled system, 
will therefore be related by 

D U) =f -2a DR (3) 

In the case that a^O, the presence of the factor f~ 2a 
implies that, although the two systems are equivalent 
from a geometrical point of view, this is not the case for 
their dynamics, i.e., for the time scales for their relax- 
ation. In order to obtain the value of the exponent a 
we proceed as follows: A characteristic time scale can 
be defined by r = tp/u>, where tp is a typical angular 
displacement between consecutive collisions and u) is the 
angular velocity in the absence of obstacles. Note that, 
due to its geometric nature, ip does not depend on L 2 p 
if the finite extension of the obstacles and the width of 
the rod can be neglected. In contrast to this, the an- 
gular velocity is not invariant under scaling, as we show 
now. If we approximate the discrete string of beads by a 
continuous rod whose mass is distributed uniformly with 
density A = M/L, one obtains for the moment of inertia 
of the rod I = XL 3 / 12. From the equipartition theo- 
rem for the energy we then find that (uj 2 ) = Ylk^T / XL 3 , 
i.e., uj oc L~ 3 I 2 . Therefore, the time scales for the orig- 
inal system and that of the rescaled system are related 
as t^) = / 3 ^ 2 t, and thus we find for the corresponding 
rotational diffusion constants = f~ 3 / 2 D R . From 

Eq. (3) we obtain therefore a = 3/4, and Eq. (2) is hence 
transformed into the scaling law: 

AtP" 3/4 = k R (L 2 p)-P , (4) 

where the dynamic exponent and the prcfactor k R are 
the only non-trivial quantities. The reliability of the scal- 
ing law (4) is controlled by two factors. One the one 
hand, it depends on the accuracy of the factorization 
Ansatz of Eq. (2). Higher order logarithmic or power 
corrections can be present when deducing diffusion con- 
stants from reptation ideas [16-20] or, from first-principle 
theories, as has been shown for the translational diffusion 
constant in the original Lorentz gas in the framework of 
kinetic or Mode-Coupling-Theory (MCT) [2,4]. In the 
probable case where such corrections are present, the ex- 
ponent must be seen as an effective exponent intro- 
duced as a fit parameter to describe a more complicated 



functional form of D R (L,p). On the other hand, it also 
depends on the accuracy of the approximation that the 
diameter of the obstacles and the width of the rod can be 
neglected. One certainly must expect that this approxi- 
mation is no longer valid once some characteristic length 
scale of the host medium, such as the typical distance 
between neighboring obstacles, becomes comparable to 
the thickness of the rod and obstacles. 

Fig. 3 shows the calculated values of D R p~ 3 / 4 for the 
different investigated lengths and densities as a function 
of L 2 p. It can be seen that, within a factor 300 in density, 
and for all the rod lengths whose relaxation times can 
be accessed in the time window of the simulation, the 
scaling law (4) is nicely fulfilled. Furthermore, we see 
that the obtained master curve shows a marked crossover 
at L 2 p = (i « 1.4, from an exponent = 1.8 to = 
2.7, thus indicating a change in the rotational relaxation 
dynamics when the length of the rod becomes larger than 
{Qi/p) 1 / 2 . In the next section we will give a physical 
interpretation of this dynamic crossover. 

Similar arguments can be used to obtain the L— and 
p— dependence of the translational diffusion constant of 
the center of mass. For this we introduce an Ansatz that 
is analogous to Eq. (2): 

Am - k CM pV{L 2 P f . (5) 

The same arguments that led to Eq. (3) can now be used 
to obtain the relation between the translational diffusion 
constant of a given system and the one for a system which 
has been obtained by rescaling the former by a factor /: 

€. = r 27 ^M . (6) 

The corresponding characteristic time scale is defined as 
t = l/v, with I a typical distance between consecutive 
collisions and v the velocity of the center of mass of the 
rod in the absence of obstacles. From the equipartition 
theorem for the energy we find (v 2 ) = 2fceT/A£, i.e., 
v <~ L~ 1/>2 . Therefore, the characteristic time will scale 
as T W = f 3/2 T, leading to A~m = / 1/2 Am- Thus, from 
Eq. (6) we have 7 = —1/4, and Eq. (5) is transformed 
into the scaling law: 

AmP 1/4 = k C M(L 2 p)- s . (7) 

Fig. 4 shows the calculated values of DqmP 1 ^ 4 for the 
investigated densities and rod lengths (left set of data). 
As for Dr, a crossover is observed at L 2 p = (\ = 1.4, in- 
dicating a change also in the translational dynamics for 
rods longer than (d/p) 1 / 2 . Though a systematic shift is 
obtained for increasing p, a description in terms of the 
scaling law (7) is reasonably good for p < 10~ 2 , where 
values of Am/9 1 ^ 4 change by less than a factor of 2 for 
a variation of density by more than a factor of 100. If 
one considers only the range L 2 p < £1, also the data 
for p = 3 • 10~ 2 follows the master curve quite nicely. In 
any case, it seems that the validity of the approximations 
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used to obtain Eq. (7) -factorization Ansatz and neglect- 
ing thickness of rods and obstacles- is more limited for 
the translational dynamics than for rotations. 

We also mention that a better collapse of the data is 
obtained if we replace p 1 / 4 by p 35 , see right set of data in 
Fig. 4. However, we do not really give an important phys- 
ical meaning to this latter value of the exponent since it 
could just be the result of taking into account a correc- 
tion to the trivial value 1/4 by lumping together finite 
size effects or higher order contributions to the right side 
of Eq. (7). 

Except for the highest investigated density p = 3T0~ 2 , 
the data in the range of L 2 p from £i = 1.4 to (2 = 22 
also seem to follow the scaling law (7), as can be seen in 
the inset in Fig. 4, where the dynamic exponent 6 = 0.4 
is obtained from the best fit that fulfills continuity with 
the law 5 = 0.7 obtained for L 2 p < (1. Finally, as shown 
by the breakdown of the scaling, it is clear that, at least 
for p > 6 • 10 -3 , the factorization Ansatz and/or neglect- 
ing finite size effects are no more good approximations 
for L 2 p > (2 = 22, suggesting a further change in the 
translational dynamics for rods longer than (^/p) 1 ^ 2 - In 
the next section we will come back to this observation. 

Some insight about the coupling between the rotational 
and translational degrees of freedom can be obtained 
from investigating the validity of the Stokes-Einstein- 
Debye relation Dcmtr = c, where c is a constant and tr 
the rotational relaxation time [21]. We define this lat- 
ter quantity as the time where the correlation function 
(cosA<j>(t)} decays to e _1 . Here A(j)(t) is the angular dis- 
placement at time t from the initial orientation of the rod 
at t = 0, and brackets denote the ensemble average. First 
we deduce a scaling law for tr, by introducing a factor- 
ization Ansatz tr = fc Tli p a (L 2 p) v . Thus, it is found that 
= f~ 2a TB,, where is the rotational relaxation 
time in the system that has been scaled by a factor of /. 
As argued above, we also have the relation t r ^ = / 3 / 2 tr. 
Therefore a — —3/4 and we obtain the scaling law 

r R p 3 /^k TR (L 2 P r • (8) 

As can be seen in Fig. 5, data fulfill the latter equation, 
with a crossover from an exponent ?7 = 0.7tof7 = 2.7 
at L 2 p = (1 = 1.4, i.e., at the same value of L 2 p as the 
crossover observed for the rotational diffusion constant, 
again indicating a change in the rotational dynamics for 
rods longer than (Ci/p) 1 ^ 2 - 

From Eqs. (7) and (8) we now obtain a scaling law for 
the product -Dcmtr: 

DcMT R p = k SE v(L 2 p)*>- s . (9) 

Thus, if for a given density p the Stokes-Einstein-Debye 
relation is fulfilled, the dynamic exponents S and i], for 
the scaling laws of respectively Dcm and tr, must cancel 
each other. This cancellation is indeed fulfilled within 
the error bar for L 2 p < (1 as shown in Fig. 6. (Recall 
that we have obtained S = 0.7 from Fig. 4, and 77 = 0.7 



from Fig. 5.) Thus, the Stokes-Einstein-Debye relation is 
violated for rods longer than (£i/p) 1//2 , in analogy with 
the phenomenon of enhanced translational diffusion ob- 
served in supercooled liquids -with temperature as con- 
trol parameter-, close to the glass transition [22,23]. 

V. DISCUSSION 

From the investigation of the rotational diffusion con- 
stant and relaxation time, we have observed two differ- 
ent dynamic regimes for rods shorter and longer than 
a crossover length l\ = (d/p) 1 ^ 2 , with (1 = 1.4. i.e., 
i\ = \.2p~ x l 2 . This crossover has also been observed 
for the translational dynamics. Moreover, for the latter, 
a second crossover, apparently not present for the rota- 
tional dynamics (see Fig. 3), has been obtained. The 
corresponding crossover length is £ 2 = (( 2 /p) 1/2 , with 
C 2 = 22, i.e., l 2 = 4.7/j- 1 / 2 . For the density p = 10~ 2 , 
corresponding to the trajectories represented in Fig. 1, 
the crossover lengths are l\ — 12 and £ 2 = 47. Therefore, 
the latter trajectories correspond to rods in the three dy- 
namic regimes separated by i\ and £2- 

In order to understand the physical meaning of the 
observed crossover lengths, we try to relate them with 
some characteristic length present in the host medium. 
For this, it must be taken into account that a set of points 
distributed randomly in the plane does not show a homo- 
geneous configuration, as e.g., in a liquid. On the con- 
trary, it consists of clusters of close points and big holes, 
see Fig. 7. It can be expected that the size and shape 
of these holes are related to the relaxation dynamics of 
the rods at short and intermediate time scales since they 
will determine how far the rod can propagate without 
hitting an obstacle. Typically such holes are not circular 
but rather elliptic, as can be recognized from Fig. 7. Due 
to the strong disorder present in the configuration of the 
obstacles, one should expect a wide distribution for the 
size of the holes and in the following we present a well 
defined method to determine such a distribution. The 
algorithm goes as follows: 

1) Pick a random point Q in the plane. (Note that with 
probability one this will not be an obstacle.) 

2) Look for the three nearest obstacles Oi, 02,03 that 
enclose Q. 

3) Select (randomly) one of the three sides OiOj of the 
triangle O1O2O3. 

4) Choose randomly one of the remaining obstacles, 
O ncw , with the conditions that i) The triangle OiOjO nev/ 
does not enclose any other obstacle and ii) The polygon 
Oi02 030 ncw has only internal angles ip < tt. 

5) If i) and ii) are fulfilled, the initial triangle is expanded 
to the new polygon O1O2O3O4, with O4 = O ncw . If not, 
O ncw is removed as a candidate to expand the initial tri- 
angle from the side OiOj, and we repeat step 3) and 4) 
until we are able to expand the initial triangle. 
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We now repeat steps 3) — > 5) (note that now in step 
3) we have a polygon of order n > 3 instead of a tri- 
angle) until i) and ii) are not fulfilled by anyone of the 
remaining obstacles and the polygon 0\02---O n cannot 
be expanded any more. Thus this procedure has gener- 
ated the "largest" possible convex polygon that has Q in 
its interior and obstacles at its corners. Note that each 
point Q will have several such polygons, since steps 3) 
and 4) contain a random element. 

To characterize the shape of each of these closed poly- 
gons we calculate its tensor of inertia (by taking mass 
unity for each corner of the polygon). The two eigen- 
values of this tensor can be used to define the ellipse of 
inertia, which in turn can be used (by means of ^4. max /2 
and A m i n /2, its major and minor semi-axis, respectively) 
to characterize the shape of the elliptic hole associated 
with the polygon. 

An illustration of this procedure is given in Fig. 7. 
As mentioned above, the resulting polygon (and asso- 
ciated ellipse of inertia) to a given point Q are generally 
not unique, and different random selections of the lower- 
order polygon sides and of the remaining obstacles in the 
consecutive upgrading steps can lead to a different final 
polygon. This arbitrariness finds its corresponding part 
in the dynamics of the rods in that there the trajectory 
depends not only on the initial location of a rod but also 
on the direction of its initial velocity. However it must be 
stressed that we are investigating quantities as diffusion 
constants or relaxation times that are obtained as aver- 
ages over an ensemble of rods. Thus we will try below 
to relate the observed features for these latter quanti- 
ties with the average size of the elliptic holes in the host 
medium. 

To measure the size and shape of the local el- 
liptic holes we used, for a large number of random 
points Q, and a few times for each Q, the above de- 
scribed algorithm. This allowed us to obtain the dis- 
tribution of the minor, g m in(A m i n ), and major axis, 
<?max(^4max)) which are shown in Fig. 8 for a den- 
sity p = 1CT 2 . The distribution g m i n (A m i n ) is rather 
asymmetric and shows a noticeable tail. It can be 
well described (see Fig. 8) by a Gamma distribu- 
tion g m in(A min ) =_(2/e) z (A min / J 4 min ) z " 1 exp[-z(A min - 
^minVAninKr^Amin)- 1 . In this expression A min is the 
average minor axis, z is a dimensionless shape parameter, 
and r(z) is the Gamma function for z [24]. The fitting 
procedure yields the values A m ; n = 12 and z = 2.7. The 
standard deviation, given by a = A m i n /^/z, is a = 7. 
On the other hand, the distribution of major axis is 
rather symmetric and can be well described by a Gaus- 
sian g max (A max ) = exp[-(A max - A max ) 2 /2cr 2 ]/\/27ro', 
with an average value A max = 47 and a standard devia- 
tion a = 14. Obviously, for any other density p the corre- 
sponding average axis -and analogously for the standard 
deviations- will be simply obtained by a scaling relation 

A(p) = Aoipo/p) 1 / 2 (10) 

where A corresponds to the -minor or major- average 



axis for the reference density po = 10 2 . Thus from 
Eq. (10) we find the relations A m i n = 1.2p -1 / 2 and 
A - 4 7/?- 1 / 2 

We have shown that there are two characteristic length 
scales in a random host medium: the averages of the mi- 
nor and of the major axes of the elliptic holes. These 
latter values nicely match respectively the first, £\, and 
second crossover length, £ 2 , observed for the diffusion 
constants. From these correspondences we obtain imme- 
diately a simple physical picture for the dynamics of the 
rod. The system has one rotational degree of freedom and 
two translational ones. The latters can be decomposed 
into components parallel to the minor and major axes 
of the local elliptic hole. For rods of length L < A m [ n , 
rotations and translations are, on average, hindered only 
weakly by infrequent collisions with the obstacles. There- 
fore, the trajectories will look like those of a typical Brow- 
nian particle, as has been shown in Fig. 1 for L = 3, and 
the translational and rotational dynamics will be basi- 
cally isotropic. 

When crossing the first characteristic length A m ; n , the 
rod, on average, will not be able to pass transversally be- 
tween obstacles that have a distance less than ^4 m i n , and 
many of the relaxation channels that correspond to a mo- 
tion perpendicular to the minor axis will be suppressed 
due to collisions. At the same time also large angle rota- 
tions will be strongly suppressed. These effects will result 
in the observed dynamic crossover seen in Dcm and Dr. 

As can be recognized in the trajectory for L = 39 in 
Fig. 1, no big blobs are present, in contrast to that for 
L = 3. Certainly, due to the hindrance on the transversal 
motions, rods of the former size cannot perform narrow 
closed loops, and big blobs are not formed in the time 
scale of the simulation. The breakdown of the isotropy 
of the dynamics also explains why the Stokes-Einstcin- 
Debye relation is violated for sufficiently long rods. For 
supercooled liquids close to the glass transition temper- 
ature, such a violation is explained as the result of the 
formation of fluidized domain bottlenecks which are effi- 
cient for the translational diffusion of the molecules but 
not for rotational relaxation [25]. A similar picture can 
be established for the generalized Lorentz gas: though 
transversal motion is strongly hindered, the rod can es- 
cape from the local cage formed by the neighbors by a 
longitudinal motion, -i.e., by moving through a bottle- 
neck between neighboring holes. Though this mechanism 
will result in a significant spatial decorrelation, it will not 
change very much the orientation of the rod, and a com- 
plete angular decorrelation will take place only for a long 
displacement from the initial position. As a consequence, 
translational dynamics will be strongly enhanced in com- 
parison with rotations, leading to the breakdown of the 
Stokes-Einstein-Debye relation. 

As there is only one rotational degree of freedom, no 
further crossovers will be observed for Dr before the ex- 
pected glass transition at large L 2 p. However, a second 
crossover will be observed for Dqm when crossing A m&x , 
due to suppression of relaxation channels through the 
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major axis. Since now rods are on average longer than 
the holes, transversal motion will be strongly hindered, 
leading to the quasi-linear paths observed for L = 79 
in Fig. 1. It is worthwhile to remark that within this 
physical interpretation, in the general n-dimensional case 
where holes are hyperellipsoids, n dynamic crossovers 
should be observed for Dqm -though it would probably 
be difficult to separate them except for very small den- 
sities. In contrast to this the n — 1 rotational degrees 
of freedom (around the n — 1 principal axis perpendic- 
ular to the rod) are equivalent by symmetry and hence 
only a single crossover should be observed for also in 
the general n-dimcnsional case. Experimental and com- 
putational tests of these predictions would certainly be 
interesting. 

Finally, as we have mentioned above, for sufficiently 
large values of L 2 p, the system should certainly show a 
"glass transition" , -in the sense of a transition from an 
crgodic to a non-ergodic state of zero diffusivity- both 
for the translational and for the orientational degrees of 
freedom, as it is observed for the formers in the original 
Lorentz gas. It is not clear whether the observed break- 
down of the scaling for the translational diffusion con- 
stant at the largest values of L 2 p corresponds to the be- 
ginning of a decay to zero through some MCT-like power 
law Dqm ~ (Cc — L 2 p) x [4]. On the other hand, for the 
rotational diffusion constant no significant breakdown of 
the scaling has been observed, suggesting that for rota- 
tions the L 2 p-investigated range is still far from the glass 
transition. Simulations at larger values of L 2 p should 
shed light about this questions. Work in this direction is 
in progress. 

It must be remembered that the results presented in 
this paper correspond to a random configuration of disk- 
like obstacles that can overlap. If we define a charac- 
teristic distance d = p~ 1//2 , we find that d is noticeably 
larger than the soft-disk diameter: d > 5a. Therefore, 
for the investigated densities, we do not expect signif- 
icantly different results if the configuration of the host 
medium was slightly modified such that the obstacles do 
no longer overlap (e.g. by displacing them by a small 
amount). (Differences must be expected at significantly 
higher densities, since for the Lorentz gas of hard disks, 
N = 1 , the details of the occuring glass transition do in- 
deed depend on the nature of disorder [26].) However, if 
the nature of the host medium is modified in such a way 
that the location of the obstacles corresponds to a liquid- 
like configuration, i.e. it shows a new intrinsic length 
scale, it might be expected that the resulting dynamics 
differs from the one discussed in the present paper. Work 
in this direction is in progress. 

VI. CONCLUSIONS 

We have investigated, by means of molecular dynam- 
ics simulations, the dynamics of a rod in a random frozen 



environment for a wide range of rod length and density of 
obstacles. The calculated diffusion constants for rotation 
and translation of the center of mass can be scaled onto 
a master curve that has been deduced from simple scal- 
ing and kinematic arguments. The abrupt crossovers ob- 
served in the dynamic exponents take place when the rod 
length crosses two characteristic length scales: the aver- 
age minor and major axis of the elliptic holes present in 
the random host medium. In analogy to the observation 
for supercooled liquids close to the glass transition, the 
Stokes-Einstein-Debye relation is violated for rods longer 
than the average minor axis of the elliptic holes. It is 
found that this effect can be understood by the longitu- 
dinal motion of the rod that leads to a fast decorrelation 
of the translational degrees of freedom but is inefficient 
for the rotational one. 
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FIG. 1. Typical trajectories of the center of mass of a rod 
for L = 3, 39 and 79 at the density p = 10~ 2 . All them cor- 
respond to a simulation time i s im = 10 6 , with a time interval 
of 500 between consecutive plotted points. The arrow at the 
bottom indicates the length scale. 
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FIG. 2. Typical trajectories of the center of mass of the 
rod for L = 7, p = 3 • 10~ 2 and for L = 199, p = 10~ 3 . 
The simulation time and time interval between consecutive 
plotted points are the same as in Fig. 1. 
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FIG. 3. Scaling law for the rotational diffusion constant. 
The solid lines are fits to power laws. The dynamic crossover 
at L 2 p = £i is indicated. 
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FIG. 4. As Fig. 3 for the translational diffusion constant o 
the center of mass. A double t/-axis representation is used fo: 
clarity. The inset shows a magnification of the data DcMp 0,35 




FIG. 5. Scaling law for the rotational relaxation time. 
Solid lines are fits to power laws. 




FIG. 7. An illustration of the method of construction of 
the elliptic holes in the random configuration of the obsta- 
cles (see text for details). The numbers indicate the order in 
which each triangle is added until the final polygon, defined 
by the thick segments, is formed. The corresponding ellipses 
of inertia are also shown. 
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FIG. 6. Scaling law for the product DqmtrP and viola- 
tion of the Stokes-Einstein-Debye relation at large L 2 p. The 
dashed horizontal line marks a constant value. The solid line 
is a fit to a power law. 
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FIG. 8. Points: Distribution of the minor and major axis 
of the elliptic holes for a density p = 10~ 2 . Lines are fits to a 
Gamma and to a Gaussian function, respectively. 
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